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Abstract 

We derive a formula for the level spacing probability distribution in 
quantum graphs. We apply it to simple examples and we discuss its 
relation with previous work and its possible application in more general 
cases. Moreover, we derive an exact and explicit formula for the level 
spacing distribution of integrable quantum graphs. 
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1 Introduction 

One of the major discoveries in the field of quantum chaology is the existence of 
universal statistical fluctuations in the spectrum of systems that are classically 
chaotic in the limit h 0. These statistics are well described by random matrix 
theory (RMT) in which the Hamiltonian of the specific system under consid- 
eration is replaced by a matrix where each element is an independent random 
variable except for global symmetries required by the Hamiltonian ||l[. Beside 
the universal aspects, some statistical properties may also depend on the partic- 
ular system under consideration. The main tool to study all these phenomena 
is the Gutzwiller trace formula that gives a semiclassical approximation to the 
density of states in terms of the periodic orbits of the corresponding classical 
system The application of this formula has satisfactorily explained some 
statistical properties that agree with RMT for chaotic systems Neverthe- 
less, there is no satisfactory complete explanation yet for the universal random 
character of the spectrum appearing from a specific Hamiltonian. 

Recently, Kottos and Smilansky have studied very simple quantum systems 
called quantum graphs that display statistical spectral fluctuations belonging to 
the class of systems with a chaotic classical limit ^, |^ . A remarkable aspect of 
the quantum graphs is that there exists an exact trace formula that expresses 
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the density of states in terms of the periodic orbits of the corresponding classi- 
cal dynamics in a similar way as the Gutzwiller formula does for Hamiltonian 
systems. These nontrivial features of these extremely simple systems have made 
them natural toy models of quantum chaology. 

In the same perspective, several papers have been very recently devoted to 
these systems |6[ ^, ^ . On the one hand, Kottos and Smilansky have studied 
scattering processes in quantum graphs showing that these systems display all 
the features which characterize quantum chaotic scattering ||^. On the other 
hand, the analysis of the statistical spectral fluctuations on graphs has been 
considered by Schanz and Smilansky |^ as well as by Berkolaiko and Keating 
[H- These last authors have studied for star graphs the two-point correlation 
function, a quantity which reflects the long-range spectral correlations. In their 
analysis they introduce ensemble averages (for example over the lengths of the 
bonds) in order to get a formula which is exploited by a combinatorial analysis. 

In the present article, our aim is different in two main aspects. Firstly, 
we want to consider the spacing probability distribution, which reflects short- 
range spectral correlations and, secondly, we want to study the dependence of 
this distribution on the parameters of the system, in particular, on the bond 
lengths. Accordingly, we do not introduce external average but we develop a 
method based on ergodicity. With this purposes, we derive a general formula 
for the level spacing probability distribution in quantum graphs using a very 
simple ergodic theorem. This formula applies more generally, too every system 
with levels determined by the zeros of a quasi-periodic secular equation. The 
result being exact, it contains all the information on the particular system. To 
obtain the universal behavior observed in some graphs from this result, further 
assumptions and simplifications should be made. We do not address here this 
difficult problem. Instead we apply our result to very simple graphs, which 
nevertheless gives interesting results (such as level repulsion) and which can 
guide the approach to more difficult and interesting cases. 

In Section^, we review some results about the quantum mechanics on graphs. 
In Section ^, we derive our main result, which is a general formula for the level 
spacing probability distribution given in terms of a Poincare mapping defined in 
a certain surface of section S. In Section U, we use the density of states for graphs 
to obtain information about S. In Section]^, we illustrate our result with some 
simple graphs. In Section ^, we compare the level spacing distribution obtained 
numerically for a complex graph, with the result of RMT. Then, in Section ^ 
we compare our result with a related theory proposed by Berry. Conclusions 
are drawn in Section ^. 

2 Energy levels of quantum graphs 

In this section, we introduce the main results known about the energy levels of 
quantum graphs in order to be complete. We refer to the works of Kottos and 
Smilansky for details[|j. 

Graphs are vertices connected by bonds. Each bond b = connects two 
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vertices, i and j. On each bond b, the component ^f, of the total wave function 
\I/ is a solution of the one-dimensional Schrodinger equation. Here we consider 
the time reversible case {i.e. without magnetic field) 

where k is the wavenumber. Moreover, the wave function must satisfy boundary 
conditions at the vertices of each bond (i and j in the previous equation) , which 
ensures continuity and current conservation, i.e., 

*fc(0) = 

for all the bonds b which start at the vertex i and 



for all the bonds b which end in the vertex j. The length of the bond b is denoted 
by lb or The current conservation reads 



= Ai ipi 



where denotes a summation over all the directed bonds which have their 
origin at the vertex i. These conditions guarantee that the resulting Schrodinger 
operator is self-adjoint. Note that, in this formulation, each bond has two 
directions and we have to distinguish between the two different directions of a 
bond. This means that the dimension of the vector ^E" = [\E'i(a;), . . . , ^25(2^)]"^ 
is 2B where B is the number of bonds of the graph. 

When Ai — > 00 (Dirichlet boundary conditions) the graph becomes a union of 
noninteracting bonds. These are called "integrable graphs" because the classical 
dynamics corresponds to particles bouncing in the bonds leading to a phase 
space with the topology of a torus. We come back to this case in Subsection 



5.4, For finite Ai, the asymptotic properties of the spectrum become independent 
of A at high wavenumbers and, indeed, there is a convergence to the Neumann 
limit where all the A's are equal to zero. In what follows we consider this case. 

As a result of the boundary conditions, we get the secular equation which 
can be written in the following equivalent ways 

det[/-S'(fc)] =0 (1) 

with S = D{k)T a unitary matrix of dimension 2B where 

Dab = Sab e^'"'" , with la = h (2) 

and 

Tab = -S^b + ^ (3) 
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if the bonds a and b are connected through a vertex (here called i) and zero 
otherwise. The notation b defines the reverted b bond. 
The secular equation can also be written as 

det h{k) = 

where h is a matrix of dimension V {V is the number of vertices in the graph) 
given by 

\ (sinfc;(,,,))-ia, if 

Cij being the connectivity matrix with elements equal to one if the vertex i is 
connected to j and zero otherwise. 

It is clear from both secular equations that the eigenvalues are given by the 
zeros of an almost-periodic function. 

Using Eqs. (|l]), and (||), it is possible to write the quantization condition 
in terms of the zeta function 



= 



p 



where p denotes a periodic orbit, Up is its period, Lp is its length, 7p is related 
to the stability of the orbit and fip is the analogue of the Maslov index. Note 
that Lp = rriili where the mi are integer numbers. If we define Xi = kli we 
can see that C,{k) = C,{xi = kli, . . . ,xb = kls) with 

C(xi, . . . , a;s) = n " e<S' '"■^'+'^''^)] (5) 

p 

Note that (^{xi,... ,xb) is 27r-periodic in each of the variables, so that C(fc) 
is an almost-periodic function. It can happen that the lengths of the graph 
are not all incommensurate. In that case, it is convenient to define a new 
function F(xi, . . . , Xn) where n is the number of incommensurate lengths, which 
gives C(fc) when evaluated in xi ~ kli, ... ,Xn = kin (here h, . . . ,ln are the 
incommensurate lengths) , i.e., 

F{xi ^kli,... ,Xn = kin) = f{k) = C{k) 



3 Level spacing distribution for almost-periodic 
functions 

3.1 The level spacings as the first-return times of a Poincare 
mapping 

In this section, we derive the probability distribution for the spacing between 
the successive zeros of an almost-periodic function /(fc). Let us call {fci};^o 
ordered solutions of f{k) = 0. 
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The probability of having two successive zeros at a distance (s, s + ds) is 
given by 

, , #{ki <K:s< ki+i -ki<s + ds} 
F(s)ds = Inn ; 

or equivalently by 

Pis) = lim ^ E '^[^ - (^'+1 - ^')] (6) 

1=0 

By the definition of /(fc), there exists a function F{xi,X2, ...,a;„) such that 

f{k) = F{xi = kli,X2 = kl2, Xn = kin) 

where the parameters Zi,^2, •■•j^n are incommensurate real numbers, which, for 
the case of graphs, form the set of incommensurate lengths, and from which all 
the other lengths can be obtained by linear combinations with rational coeffi- 
cients. 

The function F(xi,X2, ■■■,Xn) is periodic in each of its arguments Xi with 
a prime period Pi. Accordingly, we can consider the function _F on a torus 
T" : < a;, < with i = 1, ...n. 

The equation 

F(xi,X2, ...,x„) = (7) 

defines a hypersurface S on T". 
The equations 

^=h = ..,n) (8) 

define a fiow on this torus. In Eq. (||), the wavenumber k plays the role of 
the time. Because of the incommensurability of the "frequencies" U this fiow 
has the remarkable property of being ergodic. We will exploit this property of 
dynamical systems theory to obtain the desired expression for the level spacing 
probability distribution. 

First, we note that each intersection of the trajectory {xi = fc^ij^^^with the 
surface E gives a zero kj S {^ili^o- Therefore, this surface plays the role of a 
Poincare surface of section for the present dynamical system. 

In this hypersurface of section, the fiow induces a Poincare map 

kn+l = fcn + T(Cn) 

where ^„ is a point on S that is mapped by the flow on fn+i also in E. These 
two points are the successive intersections of the trajectory {xi = kli}^^^ with 
the surface E at the times fc„ and fc„+i respectively, t(^) is the time of first 
return to the surface of section E. 
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From Eqs. we have that 

kn+l - kn = t[?!)"(^o)] 

in which is an initial condition belonging to S where the iteration started. 
Now, we can write the spacing probability distribution (|^) in the form 

1 

P(.^) - nT.H'- r[^'i^o)]} (10) 

The ergodicity implies that the value of the distribution ( p^ ) is almost every- 
where independent of the initial condition so that can be any point on 
the torus T" and not necessary one corresponding to a zero. Moreover, the 
ergodicity implies the existence of a measure z/ on E which gives the spacing 
probability distribution according to 

P{s) = / S[s - t(0] (11) 

We now turn to the determination of this invariant measure v. 



3.2 The invariant measure u 

When the lengths h are incommensurate, the dynamical system (^) is ergodic 
on the torus. That is: For any measurable function g{xi, ...,a;„) defined on the 
torus we have that 

1 '■^ 



lim - / g[ip'{xo)]dt^ ^i{dx)g{x) (12) 

where (p*{xq) = It + xq is the flow (i^*, I and xq are n-dimensional vectors) and 
fj,{dx) = is the Lebesgue measure on the torus. 

Let us define the function At[(p*{xo)] as the time of flight of the trajectory 
after the last crossing of the surface of section S. If the last crossing happened 
at kn then At[(^*(xo)] = t — kn- 

We replace the function g by 

g[ip\^o)]^Q{s-At[ip\^o)]} (13) 

and we compute in this case the integral of the left-hand side of Eq. (|l 

r g[^'mdt = Y. r Q{s~At[^*{^o)]} S{t-kn)dt 
Jq Jo 

We assume that there are N zeros in the interval [0, T] and we call them 
ko, /cjv-i so that we get 

„T N-l N-l 

/ g[^*(eo)] dt = Y, e[s - (fc„+i - kn)] = E Q - r[rm} 

•^0 Ti=0 r!=0 



6 



For large values of T we can consider that T — kpf with N the number of zeros. 
Denoting by (d) the mean density of zeros, we have N — (d) and T = so 
that we finally get 



hm i / g[vp*(eo)] dt = id) hm ^ ^ ^ " (Co)]} 



(14) 



We recognize the cumulative function times the mean density in the right-hand 
side of the expression (p^). 

To compute the right-hand side of Eq. (|l^) , we have to write g as a function 
of the coordinates x. For this purpose, we remember that 



^5(t-fc„) = |/'WI^[/(t)] 
{»} 



(15) 



Now j[t) — F[y*(go)] and fit) = VF[(^*(^o)] • I. Replacing these expressions in 
Tsh, and ([151) in (flSh, we obtain 



<?[v'*(eo)] = e{s-A%*(eo)]} \^F{ip\i^)\-l\ 5{F[^*(eo)]} 

from which we infer 

g{x) ^ Q\s ~ l\t{x)\ |VF(a;)-;| &{F{x)\ 

where ts.t(x) is the time taken by the trajectory to arrive at x since its last 
crossing with S. 

Now, we compute the right-hand side of Eq. (|l2|) which we denote by /: 



liidx) gix) 



1 



IT" 



dx g(x) 



(16) 



We perform the nonlinear change of coordinates a; — > (i, i^) where ^ are the n — 1 
coordinates that parametrize the surface S, i.e. 



Xi = lit + s.i{i) 



(17) 



where the functions Si(^) are such that F[si(^), s„(^)] = 0. In the new 
coordinates, the equation for the surface E is i = or t = t(^). In these new 
coordinates, we have that 



At(x) = t 

dx = J{i) di dt 



with the Jacobian determinant 



h ■ 


In 


dsi 








dsi 









(18) 
(19) 



(20) 



7 



and < i < r(^) where t(^) is the time of first return previously introduced. 
In these coordinates, / is given by 

^^W\IJ^ ^(6 f^^^'^ dt Q{s -t)\VF-l\ S[F{i, t)] 
The integration over t can be carried out using a new variable u defined through 

u{t)^F{^,t) (21) 

where ^ is kept constant. Differentiating with respect to t gives ^ = VF • I and 
we get 

I^T^ I d£, JiO I , „ e[s - t{u)] \VF ■ l\ 5{u) 

This integral picks up the value of t{u) at w = 0. From Eq. (|2l]), we see that 
u = is the equation that defines E and, as we noticed after Eq. (p7[), there 
are two solutions t{0) = or t(0) = r(^) in the new coordinates. But since t is 
the "time of flight" after the last crossing, we consider the second solution and 
we finally get 

I=^J^d^mQ[s~rm. (22) 

From Eqs. (p^), (p^, (p6|), and (p^), we find the cumulative function and by 
differentiation with respect to s we obtain the level spacing probability density 

N-i ^ r 

^^'^ = iToo N ^ - "['^"(^"^l^ = L '^^ ^^^^ ' "^^^^ ^'^^ 

On the other hand, the density can also be expressed in a geometrical form. 
Indeed, starting from its definition 

id) ^ lim Mhl^ = lini U'^Y. - fc") 

•^0 {/} 

and using Eq. (|l^) and Eq. we have 

{d)^^J^^^dx \VF-l\ 6[Fix)] 

Rewriting this expression in terms of the new coordinates ( p^ ) and then doing 
the changes of variables (|2^) , we obtain 

{d)^j^JjCJiO (24) 
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Let us observe that this expression (^4|) for the density can be obtained directly 
from Eq. ( |2^ ) and the normaHzation condition P(s)ds — 1. 
Accordingly, we can write the spacing probability density as 



which is the central result of this paper. The expression ( pS]) has a very simple 
geometrical interpretation. It gives the spacing probability density as the ratio 
between the flux of the probability current lS[s — t(^)] through the surface E 
and the flux of the constant velocity field I through the same surface S. 

From Eq. (|25|), wc can conclude that the invariant measure in ( pj]) is given 

by 



4 The density of states as a sum rule for graphs 

In the previous section, we derived a formula which relates the density of zeros of 
an almost-periodic function to the properties of the surface of section S defined 
in a torus. The dimension of the torus equals the number of incommensurate 
lengths and the periodicity Pi in each variable depends on the relations between 
the length li and those which are commensurable with it. For example, if there 
is a length commensurable with ^i, i.e., Ik = then the variable xi will have 
the period Pi = 2Trq. In the case where the relation is of the form Ik = nli, 
or all the lengths are incommensurable, we can always consider that Pi = 27r, 
Vi. In what follows, we consider this to be the case. As a consequence, we can 
rewrite Eq. ( p^ as 

As we have already pointed out, this expression has the geometrical interpre- 
tation of a constant flux I through the surface S. Because of the periodicity of 
E in the n-dimensional real space i?", we expect that the projection of E in all 
of the n directions covers the complete plane. (This would be false if E was a 
closed surface but we suppose that this is not the case.) Therefore, if we call Ej 
the projection of E in the i^^ direction we have 

and dsi = TOi(27r)"^^ with rrii the number of sheets of E for the projection 
in the i**^ direction. Consequently, we have 

(d) = ^^m, (27) 
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The number rrii can be determined for each particular case by inspection on the 
quantization formula [for example Eq. (|l|)] . 

Here, our purpose is to revert the argument and use this formula to obtain 
the TTii. This is possible because there is a general expression for the density. 
For graphs, the density of states was obtained in Ref. using the properties 
of (|l|), (||) and (||) together with a formula for the density of states derived in 
the approach of scattering quantization. The result is simply given by 

{d) = ^ . (28) 

TT 

From Eqs. (p7|)-(p8|), we get the desired equation for the m^'s 

J2 k = 2Ltot (29) 

i 

Since the sum is only performed over the incommensurate lengths, forming a 
basis from which all other lengths can be obtained, this equation gives all the 
rrii because we can write Ltot in such a basis. 

We can deduce from here that when all the lengths are incommensurable 
there will be two sheets in each projection. 

It is possible to reduce the "volume" of the torus by noticing that, in fact, 
we do not need the function F(xi, . . . ,Xn) to be periodic but the surface S. 
Since this surface is given by F{xi, . . . , Xn) — 0, the period Pi with which the 
surface is repeated in i?" is given by the smallest of the period or anti-period 
of the function in the variable Xi [i.e., the values Pi for which F{xi, ... ,Xi + 
Pi, . . . ,Xn) = iF(xi, ... ,Xi,... ,Xn)]. Note that we call it again Pi but there 
is no risk of confusion. Moreover, in the rest of the paper, we shall use this 
definition. 

5 Application to simple systems 

5.1 A three-bond star graph with three different bond 
lengths 

Let us consider the simple graph composed by three bonds attached to a vertex. 
The spectrum of this graph is given by the zeros of the function 

/ (fc) = cos kli cos kl2 sin kl^ + cos kli sin kh cos kl^ + sin kli cos kl2 cos kl^ (30) 

The function ( |30| ) is an almost-periodic function. Let us define 

G{xi, X2, X3) = cosa;i cosa;2 sin 0:3 -I- cos xi sina;2 C0SX3 -I- sinxi COSX2 cos 0:3 

This function is 27r-periodic in each argument but has an anti-period tt. It 
satisfies 

G{xi = kli,X2 = kl2,X3 ~ kls) = f{k) 
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The equation G{xi,X2,X3) = defines a surface S with a double cone joined 
by a singular point. The singular point is given hy xi = ^ , X2 = , X3 — ^ (see 
Fig. |l|). For simplicity, we translate the coordinate system to that point so that 
we consider the function 




Figure 1: Plot of the surface S for the three-bond star graph with three different 
lengths. The plot is obtained from G{xi, X2, X3) — 0. 

As we saw in Section ^, the shape of the surface S determines the level 
spacing probability distribution. For small spacings s, the distribution is given 
by the iterations with short "times of flight" . These are determined by intersec- 
tions near the singularity of T, because there are arbitrarily close points in its 
neighborhood. In order to study the behavior of the level spacing probability 
distribution for small spacings we carry out our analysis near the singular point 
where the function F can be approximated by the quadratic function 

F{xi,X2,X3) = X1X2 + X1X3 + X2X3 + 0{Xi) 

We diagonalizc the quadratic form with a rotation of coordinates and we finally 
get 

F(yi,y2,y3) - 2yf - yl - yj + ). (31) 
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In the y-coordinates, the flow is given by — where e\ = [li + 12 + h) / 
62 = {h — h)/V2, 63 = {I2 + I3 — 2I1)/ Vs. Now, we apply our theory. We define 
new coordinates through the transformation (j/i, 1/2, ys) (77, t) 



Vi = si(?7,$) + eit 
2/2 = S2(r/,^) + e2i 

2/3 = S3(?/,C) + 63* 



(32) 



where the functions Si{ri,^) are zeros of Eq. i.e., 2sf — S2 ~ — 0- A 

solution is 



r?2 +^2 



(33) 



Eqs. ( p2| ) and ( |33| ) define the new variables. We need to compute J and t(?7, ^). 
For J, the calculation is straightforward. Using (poh and (^), we get 



where 

Kv,0 = 2eisi 



J 



6252 - 6333 



(34) 



eiV2(r/2+e2)+ 62^-6377 +0(2) (35) 



The notation 0(2) means here "to second order in 77,^ or t" . If we write Eq. 
(^) in the new coordinates we get 



F{r,,^,t)^ah^ + 2tbi7j,0 + iH 



l) + 0{3) 



where b is defined by Eq. (|3^) and 

2 o 2 2 
a = — 62 - 



= 2(?iZ2 + Z2Z3 + ^i^3) 



(36) 



(37) 



The third term in Eq. (36) is zero by definition. In the new coordinates, the 
surface of section E is given by the roots of F{ri,£^,t) = 0, i.e., t — Q and 
t = The function t(?7,^) represents the "time of flight" of a trajectory 

which starts at one point on the lower cone with coordinates {rj, ^) and arrives 
to the upper cone. That is 



'(^'C)---^ = 



VW+F) + 62^ - 637/ 



2e2 



0(2) 



(38) 



Now we are ready to compute P{s) for small s using (|2^). As we have already 
noticed the integral in the denominator is just the density of states which is 



(d) 



h + h + h 



(39) 
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for this graph. The integral in the numerator is 

J VW+¥) V "V 

which is performed by changing to a variable u(^) = * + ^ where 77 is kept con- 
stant, by using (p8[), and then by integrating in rj. The details of this calculation 
are left to the reader. The result is 



This, together with @ and (H), gives 



'1 + '2 + '3 TT 

Usually, we express this probability density in the scaled variable A such that 
the mean level spacing is equal to unity: 

n^) - + ^(^') (40) 

(ii +12 + h) 

We observe that this simple graph already presents the Wignerian level re- 
pulsion, a property usually associated to chaotic classical dynamics. To our 
knowledge, there are only a few systems for which this result can be derived 
exactly. 

In Fig. ^, the cumulative function is depicted as a function of and 
compared with a numerical calculation of the spacing distribution. The slope 
at the origin is half of the slope of P(A). The straigth line in the figure has half 
of the slope given by (^o|). We see that there is very good agreement between 
( ^o|) and the numerical result. 

There is an interesting point about this result. The slope of ( ^o|) takes values 

between zero and ~ 0.6 as the lengths hjhJs vary. Therefore, the slope 

. . . 

only varies on a relatively small interval. This means that changing the length 

of the bonds (but always keeping them irrationally related) does not change 
very much the slope of the spacing probability density -P(A). Moreover, we 
note that the dependence on the lengths can be seen as a quotient between two 
different averages of the lengths (a geometric average and a arithmetic one). 

It is also interesting to notice in Fig. ^ that the projections of the surface E 
onto each axis cover the corresponding plane only once in the torus of volume 
TT^, which is implied by the formula ( ^9|) of Section ^ and by Eq. (^9|). 

5.2 A three-bond star graph with two different bond lengths 

Now, we consider the same graph as in the previous subsection but with only 
two different lengths, say ^i, ^2- Taking li = in (|30|), we get the function 
which gives the zeros for this graph 

f{k) = cos kli (2 sinfcZi cos + sin kl2 cos kli) 
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F 




Figure 2: Plot of the cumulative function F — P{A')dA' of the level spacing 
distribution, as a function of for the three-bond star graph. The straight 
line is the prediction obtained by integration of Eq. (|40|). Here, h = tt, I2 = 
3.183459012, and ^3 = 3.1442336073. 

We now introduce the function 

F(xi,X2) — COS ail (2sina;i C0SX2 + sinx2 cosxi) 

such that F(xi — kli,X2 — kl2) — f{k). The function F{xi^X2) is tt periodic in 
xi and TT anti-periodic in X2 and can be considered in the torus < < tt with 
i = 1, 2. In Fig. ^, we draw the lines where F{xi, X2) = in the plane (xi, X2)- 
Changing the origin of the coordinates to the singular point (f , §) corresponds 
to analyze the function 

G{xi,X2) = — sinxi (2 sin 0:12 cos xi + sinxi cosa;2) 

Note that around the singularity the function G can be approximated by the 
quadratic form 2xiX2 + xf. In this example, we explicitly obtain the density of 
states using (d) = :pT J-^"^ ^'C- First, we identify the surface E as the union of 
two lines Ei and E2 as indicated in Fig. || so that 
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Figure 3: Schematic representation of the torus for the problem of Subsection 



5.2| . The surface E is composed of the xi-axis and of the curved curve. The 
trajectories that cross the surface E belong to the region A, B, C or D. Note in 
the figure that the regions D are one next to each other that is why we put 2D. 



where Ei is given by sin.Ti = 0, i.e., Ei ~ {xi = 0, < X2 < §}• This surface 

h h 

is given by si(^) = and S2(^) — ^ so that J — dsi ds2 = ^i- The surface 

Hi ~di 

E2 is given by 2 sin 0:2 cos xi + sinxi cosa;2 = 0, i.e., X2 — — arctan tana;i). 
Therefore, if si{£) = ^ we find S2(^) = — arctan (i tan^) so that we get J = 
I 21 



l+cos2 5 



. Consequently, we obtain 



^2 + 



2/1 



1 + C0S2 ^ 



2h+l2 



as expected. In Fig. ^, we notice that there are two sheets of E with projection 
onto X2 and only one with projection onto Xi as expected from the formula ( p9[ ) 
of Section |. 

Let us compute the level spacing probability density P{s). From the symme- 
try of Fig. ^, we recognize four regions {A, B,C, D), each one repeated twice. 
For three of these regions {A,B,C), a trajectory joins a straight line with a 
curved one. For the other region (D), two curved lines are joined. In this 
respect, we need two expressions for the "time of flight": ri(^) and T2{^). 

The "surface of arrival" is determined by 2 sin 0:2 cos xi + sin cos 2:2 = 0, 
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i.e., X2 = — arctan (^tana^i). Considering X2 = ht + £, and xi = lit, we get 



^ = —l2t — arctan 



-tan(?ii) 



(41) 



Solving this equation for t = t{^), we find that ti(^) — |mini(^)|. From the 
periodicity of the arctangent function, the next solution [let us call it ^2(0] gives 
the time to cross the following surface (see Fig. ^ so that T2{^) — ^2(0 
Since the parameter ^ moves in the a;2-axis, we can write 

J 6[s - TiO] ^ 2h r S[s-niO]d^ + 2h ^ S[s - T2m 

Jo JQ 

where P = it and Q = — ^ . The first integral takes the contributions from the 
regions {A, B, C) and is easy to compute with the change of variable u = ti{S). 
We find 

du 



5[s~Ti{C}] di 



dri[£,(u)] 



6(s — u) 



where Tp 



ri(^). Differentiating HWj with respect to we get 



drii^iu)] 



d^ 



2/i 



SO that 



S[s-n{m dC 



2li 



^2 + l+3cos2 ii< 

, 



, for s < Tp 
for s > Tp 



For the second integral which takes the contribution of the region D, we have 
a similar formula but it depends on the implicit functions ti(^) and ^(s) given 
by the equation T2(^) = s. If T2{Q) = - 2ti(-^) < s < tp we get 



S[s-T2m d^^ 



h + 



2li 



l + 3cos^h[s + Ti{^{s))] 



+ <l2 + 



2h 



1 + 3cos2/i[ti(C(s))] 



and zero if s > Tp. Hence, in the scaled variable A = -^^s, we have 



^tot 



12 + 



2li 



l + 3cos2ZiAj^ 



itot 



2li 



2 r(A) 



(42) 



where r(A) = if A < ^T2{Q) and 

r(A) ^ 



12 + 



2h 



I2 + 



2h 



l+3cos2/i|ri[{(A-c^)]| 



(43) 
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if %^T2((9) < s < ^Tp. Finally, we note that P(A) = if A > -^rp. 

The fact that r(A) has an implicit dependence on A, makes difficult its 
actual evaluation. Nevertheless, in the numerical example considered below, 
the interval ^^T2{Q) < s < ^^Tp where r(A) is different from zero is small 
and we can consider a simple approximation for r(A). 




Figure 4: Numerical and theoretical level spacing probability densities for the 
graph of 3 bonds with 2 different lengths, Zi = tt and I2 = 1.53183459. 

Consider h = t: and h = 1.53183459012. Thus, we get ^T2{Q) = 1.345 and 
^=^Tp = 1.522. Therefore, r(A) is different from zero in the interval 1.345 < 



In this interval, we can consider 
we see that the dependence on ^ is 



A < 1.522 as can be observed in Fig. 
T2{^) as a linear function of ^ (from Fig 

in fact smooth). Accordingly, r(A) is simply given by the constant tj^qj 
[Remember that Q = —2^77- after Eq. (^l|).] Subtituting the numerical 



_2_ 



values, we get = 1.107 which added to the first term of Eq. (Ea) predicts a 

peak of the order of 1.8, which agrees with the peak of the numerical result shown 
in Fig. ^. A more accurate comparison can be done through the cumulative 
function F(A) = jfdA'P{A'). Fig. | shows the numerical result and the 
analytical result obtained by integration of Eq. (^) using the approximation 
YT^T = 1.107. Here, a good agreement is observed. 

Let us note that if we consider this problem but with I3 — pli, with an even 
integer p, the curve Si does not intersect S2. As a result, the level spacing 
probability density is zero between A = and a value Ac. This is ilustrated in 



17 





0.25 0.5 0.75 1 1.25 1.5 1.75 

A 

Figure 5: Cumulative function for the same graph as in Fig. ^. The sohd Hne is 
the theoretical calculation done in the text and the dashed line is the numerical 
result. 

Fig. I with p = 2. 

5.3 Other simple graphs 

Here, we study the level spacing in simple graphs with bonds connected to one 
vertex and thus forming a loop. 

First, we consider the graph formed by two loops attached to a single vertex. 
This graph has the form of an eight. The zeros are determined by the function 
/(fc) = F{xi — lik,X2 — hk) where 

F{xi,X2) — (cosa;2 — l)sina;i + (cosxi — l)sina;2 (44) 

This function is 27r-periodic in each variable and the surface E obtained by 
F(xi,X2) = is considered in the torus — tt < xi < tt and — tt < X2 < tt. It is 
easy to see that this surface is composed by Si = {xi = 0, — tt < X2 < tt}, S2 = 
{x2 = 0, — TT < xi < tt} and S12 = {x\ + X2 = 0, — vr < X2 < tt}. These three 
sheets intersect at the singular point x\ = X2 = 0. The function i^(xi,a;2) = 
can thus be replaced by the cubic form xiX2{xi + X2) = 0. The level spacing 
probability density can be written as 



18 




Figure 6: Here, we have considered ^3 = 2li in the three-bond star graph to 
illustrate the existence of a critical value at which -P(A) is different from zero. 
Here h = V2 and h = \/3. 



with Ji = h, J2 = h, J12 = I1+ h, Ti(0 = T2{0 = 77^ for -tt < ^ < 0, 
ri(e) = T^iO = ^ for < e < TT, ri2(0 = ^ for < C < ifrl^ and 
Ti2(C) = for — < ^ < 0. Performing the integrals by using the variable 
A = ^-^^s, we get 

^^^^ I ~ , otherwise ^^^^ 

In this example, the spacing probability density P(A) is independent of the 
system parameters. We have confirmed this result with numerical calculations 

(data not shown). 

Another graph of a similar type is the one composed by a bond and a loop 
attached to a vertex. This graph has the form of a nine. Here, the surface of 
section E is given by the equation 

F{xi,X2) = 2cosa;i cosa;2 — 2cosa;i — sinxi sina;2 (46) 
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The surface can be considered in the torus — 7r/2 < xi < 7r/2 and — tt < X2 < 
and it is given by 

El {a;2 = 0,-7r/2 < < 7r/2} 

S2 = {tana;i = 2(cosa;2 — l)/sinx2, — vr < a;2 < tt} 

In this example, the calculation is similar to the one for the star graph with 
three bonds of two different lengths and we do not present it here. We only 
compute P(A) in the limit A ^ 0. For the small spacings, we can consider 
the quadratic form around the singularity at (xi = 0, X2 = 0) which is given by 
F{xi,X2) — X1X2 + x\ for xi,X2 small enough. With this approximation the 
calculation is similar to the one of the previous graph. The result is P(A) 
when A 0. 

5.4 Graphs with disconnected bonds 

The formula ( p5| ) can be used to study the spacing distribution for the "inte- 
grable graphs" discussed in |^ . These graphs are obtained by imposing Dirichlct 
boundary conditions on the vertices and are called integrable because the clas- 
sical dynamics in the graph correspond to a particle that bounces in a bond in 
a periodic motion which corresponds to a torus in phase space. In this case, the 
eigenvalues are obtained by the equations 

sin fc4 = , V6 , 

i.e., 

n 

F{xi, . . .Xn) = J^sinxi =0 

1=1 

which is the equation for the surface E. This surface is composed of all the 
faces of the n-dimensional cube which defines the torus when we identify the 
corresponding boundaries, so that E — [J - E; with E^ — {xi ~ 0,0 < Xj < 

TT , «} 

In this case, the level spacing probability density ( p5| ) is given by 

= T^4r E / Jk S[s ~ Msk)] dsk (47) 

where 



Jk — I 



k 



TT - X° TT 



and Tfe(sfe) = min<^ — r^^-r) ■ (48) 
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In the Appendix, we prove that Eq. ( p7| ) together with Eq. ( p8| ) are equiv- 
alent to: 



Ik U 



e 



L 



tot 

1 



A 



h 
Ltot 



n 



(49) 



where /i is the largest length of the graph. 

The distribution ( ^9|) is in general different from the Poisson distribution. 
The Poisson distribution is the limit of (^ ) when — > 0. Indeed, in this 
limit, the delta peak vanishes, the Heaviside function equals one and, since 
Ltot, we find 



P(A) = lim TT 1 



Let us remark that this limit means that the number of bonds goes to infinity 
but the lengths are kept constant. 

In Fig. (^, we plot the distribution (49) for different numbers of bonds. We 
observe that, for two bonds, the distribution is constant (except for the delta 
peak) and that, for three bonds, it decays linearly. In Fig. (^), we plot ( p9[ ) and 
the numerical result for a graph of eight disconnected bonds. We observe the 
very nice agreement with the formula (49), as well as the convergence toward 
the Poisson distribution. 

We want to comment on the deviations with respect to the Poisson distribu- 
tion. First, we observe a maximum spacing which is easy to understand. The 
"regular" spectrum consists in a superposition of spectra (i-e., equally 

spaced levels). The largest spacing in this superposition is equal than the spac- 
ing J- where li is the largest of the lengths In the scaled variable of unit 
mean spacing, this is ^^ff^- This maximum spacing will appear repeatedly over 
the whole fc-axis creating the delta peak in the distribution. 

There is another interesting deviation with respect to the Poissonian distri- 
bution. We can see from (E9h that the probability density of finding two levels 

in coincidence is P(0) = 1 — ^ '"^'^ compute this probability 

in another way: Writing the level density p{k) in the scaled variable x of unit 
mean spacing 

■y ^ y ^ ^ / TlLtot 

j — l m—{] 



and using the Poisson formula for the Fourier transform, we obtain the "power 
spectrum" 



1 f 



EE 



du e^y- {p{x)pix + «)) = >> ^5 U - 27r-^ 



IjUl 
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0.5 1 1.5 2 2.5 3 3.5 

A 

Figure 7: Theoretical level spacing distribution for graphs with disconnected 
bonds given by Eq. (^9|). Here, P means -P(A). We have omitted the delta peak 
from the curves. The constant distribution is for the graph with 2 bonds. The 
linear distribution is for a graph with 3 bonds. Then, we plot the distributions 
for graphs with 8 bonds and with 30 bonds, respectively. The last one is close 
to the Poisson distribution that is also plotted but it starts at a smaller value 
as predicted from (|49|) . The lengths are given by the formula li — \fi except for 
= Vl67, U = \/l07, ^8 = exp(l), ^9 = Vl05, /le = \/Tl9, and ?25 = Vl34. 

where () is the average over x and p represents the fluctuations of p around 
1 (the mean density in this variable) . Now, the mean number of levels in the 
interval [x + A, 2: + A + dA] given that there is a level at x is provided by g{IS)dlS. 
with O 

/OO 
dy e^y^ 
'OO 

Thus, in the case of integrable graphs, we get 

;2 " 

i tot ^_^g 

We see that g(0) — P{0). In order to compute the level spacing distribution, 
it is often assumed that P{A) is proportional to g{A) and that the levels are 
uncorrelated, so the probability of having two neighboring levels at a distance A 
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Figure 8: Numerical and theoretical level spacing probability densities for a 
graph with eigth disconnected bonds. Here, P means P(A). We have omitted 
the delta peak from the theoretical curve (^9|) but we see that its position co- 
incides with the numerial peak. The histogram was built with 134050 spacings. 
The lengths are = 7167, h = V2, h = VS, h = VT07, h ^ k = v%, 
Ir = -s/T, and Is = exp(l). 

is given by P{ A) = (;( A)e^ g(x)dx ^ ^^^-^ ^Yig^t these assumptions 

are not justified in the case of graphs but they are approximately valid for the 
very small spacings and also for the case of graphs with infinitely many bonds 
where the distribution is the Poisson distribution. 

We have explored the dependence on the lengths of the bonds in the level 
spacing distribution (^9|). Figure ^ shows the deviations with respect to a Pois- 
son distribution for two sets of lengths and the difference between them. We 
see that the dependence on the lengths for a graph of 10 lengths is very weak in 
the integrable case. We also observe in the figure that the deviations from the 
Poisson distribution is maximum for A = 0. 



6 Level spacing in complex graphs 

We have computed the level spacing distribution for a fully connected pentagon. 
In figure |l^ we depict the cumulative function obtained numerically with more 
than 100000 levels together with the RMT prediction iQ. Although the agree- 
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Figure 9: Deviations of Eq. (^ ) with respect to the Poisson distribution for 
for a graph with ten disconnected bonds. Here, SP = exp(— A) — P{A). The 
dot-dashed line is for the set of lengths h — \/3, I2 = V5, ^3 — V7, h = '\/TT, 
^5 = VT3, le = ^/T7, h = VT9, h ^ V^, k = \/29, Iw = VsT. The soHd 
Hne for the lengths h VToT, h = \/l03, ^3 = VW, h = \/T09, ^5 = VTT3, 
h = V127, ^7 = V13T, = VT37, ^9 = Vl^, ^10 = VT49. The long dashed line 
represents the difference between the densities evaluated with (^) in the two 
different cases. 

ment is very good some systematic deviations exist. In figure |l^ we plot these 
deviations for three different sets of lengths. We can see that they are very 
close to each other showing that the fluctuations are independent of the graph 
lengths. The dot-dashed line in figure Jill represents the fluctuations around 
RMT for a fully connected tetrahedron [pp We can conclude from these results 
that the fluctuations around RMT depend on the topology of the graphs but 
does not depend much on their lengths. Moreover, we observe that the pentagon 
(a graph of 5 vertex, 10 bonds and valence 4) has bigger deviations with respect 
RMT than the tetrahedron (a graph of 4 vertex, 6 bonds and valence 3). This 
behavior is reminiscent of an observation in that for a star graph of 15 bonds 
the form factor deviates more from RMT that for a star of 5 bonds (see also 
[H ) . These results would suggest that the valence plays a role in the deviations 
with respect to RMT. 




24 



0.8 



0.6 

F 

0.4 



0.2 



- 




/ 


- 














// 








// 
















// 








/I 








// 








// 








// 








// 








// 













12 3 

A 

Figure 10: Cumulative function of the level spacing distribution for a fully 
connected pentagon. The dashed line is the numerical result for the pentagon 
with lengths Li = 0.6li and h — \/2, I2 = V^, ^3 = V5, U = V6, ^5 = V7, 
Ig = TT, ^7 = exp(l), ^8 = VlO, I9 = Vtl, ho = VT3- The sohd hue is the RMT 
result. 

7 Comparison with Berry's theory 

Berry has studied the level spacing distribution in classically chaotic systems 
with a similar idea as the one we have developed here ||l^. He noticed that, for 
a typical Hamiltonian with real eigenfunctions (which is the same situation as 
the one we consider here), it is necessary to vary two parameters in order for 
two levels to be degenerate. This is the content of a theorem originally due to 
von Neumann and Wigner. It also implies that, in the three-dimensional space 
of the two parameters A and B and of the energy E, the eigenvalue surface 
E = E±{A,B) has the form of a double cone with its sheets joined at the 
"diabolical point" {A* , B* , E*), where A*,B* are the parameters for which the 
degeneracy occurs. Following Berry, these cones are distributed in the space 
{A,B,E) according to a unknown probability distribution p{A,B,E). Berry 
has also considered a probability distribution which rules the geometry of the 
cones [7r(a, 6, c) where a, 6, c are the parameters in the quadratic form which 
defines the cone]. 

The level spacing probability distribution is given by the average over 
energy, which can be considered in the semiclassical limit where infinitely many 
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Figure 11: Deviations of the cumulative function of the spacing distribution for 
different complex graphs with respect to the RMT result: 5F = F — Frmt ■ The 
long dashed line represents this deviation for a fully conected pentagon with the 
same set of lengths as in figure |o[ The dashed line represents this deviation for 
the lengths Li = OAli with the k of the first set in figure ^ and the solid line for 
the lengths Li = 0.14/^ with the k of the second set in figure ^ The dot-dashed 
line represents this fluctuations for a tetrahedron with the lengths Li = 1.05/^ 
h = V2, h = VS, h = TT, /4 = V6, ^5 = V7, k = VTS. 



levels lie near any given E. As a consequence, the level spacing is given for small 
spacings by the successive crossings of the conical surfaces with the line {E, A = 
Aq^B = Bq) where Aq and Bq are the parameters of the actual Hamiltonian 
under study. Berry argues that, since there is nothing special about the system 
with the parameters (Aq^Bq), the energy average can be augmented by an 
ensemble average over a region {A,B) near {Ao,Bo). Whereupon, the level 
spacing becomes 

P( A) ^ Pi^o^ Bn E) r 7r(a, b, c) [ dAdB5(A~ ^/ aA^ + 2hAB + cB2) 

kAe)) J J 

After the change of variables a — A/ A, (3 — B/A, the previous equation gives 

P(A) - A 

where the proportionality factor involves a geometric average. Berry's argument 
shows that the level spacing density should vanish linearly in generic systems 
because of the level repulsion, as expected from random matrix theory. 
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The main difference between Berry's derivation and our derivation is tliat he 
introduces by hand the ensemble average. In our derivation, the ensemble aver- 
age naturally appears from a rigorous equivalence between the energy average 
and the ensemble average given by the ergodic theorem. This ensemble average 
introduced by ergodicity has the advantage of keeping all the specificities of the 
system, i.e., the dependence on the lengths of the graph. We expect that these 
specificities disappear for graphs which are sufficiently large, in a way which has 
still to be understood for graphs with connected bonds. 

8 Conclusions and discussion 

In this article, we have derived a formula for the spacing probability distribution 
of the energy levels of quantum graphs and, more generally, for systems where 
the secular equation is given by an almost-periodic function. Our formula is 
based on the ergodic properties of a continuous-time dynamical system defined 
on a torus. This ergodic flow induces a Poincare map in a certain surface of 
section which corresponds to the locus of the energy eigenvalues in the phase 
space of the flow. The level spacings are explicitly related to the times of first 
return in the surface of section. The level spacing distribution is thus given 
by the distribution of the first-return times of the ergodic flow in the Poincare 
surface of section. 

We have applied this formula to different graphs. In general, the slope of 
the spacing density -P(A) at A = depends on the system parameters and we 
have been able to calculate explicitly this dependence in several graphs. 

We have also studied in detail the "regular" spectrum of intcgrable graphs. 
One important application of our formula ( p5[ ) is the following 

Theorem: 

If the bonds of the graph are disconnected so that the spectrum is a superpo- 
sition of n independent equally spaced spectra of wavenumbers and if the bond 
lengths are mutually incommensurable, the level spacing probability distribution 
is exactly given by Eq. when the distribution is expressed in the variable 

where the level density is equal to one. The distribution converges to the 

Poisson distribution in the limit n ^ oo. 

On the other hand, for large connected graphs, the level spacing distribu- 
tion is close to the Dyson-Gaudin-Mehta spacing distribution of RMT although 
deviations are numerically observed which depend mainly on the topology of 
the quantum graph. The deviations with respect to RMT arc more important 
for smaller graphs than for larger graphs but the effect of Wigner repulsion is 
still present in very small graphs where we observe that the spacing density also 
vanishes linearly like P(A) ^ A. 

The different results we have obtained can be understood on the basis of 
the general properties of the surface of section S, which plays a particularly 
important role. First of all, we remark that the surface E is defined as the 
set of the zeros of f{k) = F{kli, ...,kln) = in the n-dimensional phase space 
{xi = kli, ...,a;„ = kin) of the ergodic flow. Therefore, the surface E is of di- 
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mension n — 1 in this space. According to the von Neumann- Wigner theorem, 
two zeros are generically degenerate only if two constraints are imposed on the 
parameters of the systems which are here the lengths li of the bonds. Conse- 
quently, the dimension of the subset of these degeneracies is n — 3, generically. 
In the following, we refer to this subset as the singular manifold. 

The aforementioned generic situation is already encountered in graphs with 
three incommensurate lengths for which the phase space of the ergodic flow is 
of dimension n = 3, the surface of section S of dimension n — 1 = 2, and the 
deg eneracy subset of dimension ti — 3 = 0. Indeed, in the example of Subsection 
5.1, the surface S forms a cone with a self-intersection at a point. This example 



shows that a spacing density vanishing linearly is generically possible as soon 
as there are three incommensurate lengths. 

However, for graphs with only two incommensurate bond lengths, only two 
behaviors are generic: either (1) P(0) ^ or (2) P(A) = for < A < Ac. For 
such graphs, the torus is two-dimensional and the surface giving the eigenvalue 
is one-dimensional, i.e., a line on the torus. Generically, this line may intersect 
itself leading to the case (1), or it may have no intersection leading to the 
case (2). Therefore, generically, we should not expect a spacing density which 
vanishes linearly like f'(A) ~ A for a graph with only two incommensurate 



lengths. This result is illustrated with the examples of Subsections 5.2 and 5.3 



In the two dimensional examples, we notice that the singular manifold has 
the dimension n — 2 and is a point on a two-dimensional torus and, as a corollary, 
P{A) starts with a finite value. This last result is supported by the fact that, 
for the non-generic integrable systems, two levels can come in degeneracy by 
varying only one parameter. Indeed, the singular manifold is of dimension n — 2 
for the disconnected "integrable graphs" because the surface S is composed by 
the faces of the cube and their intersections are of dimension n—2, as we saw in 
Subsection 5.4. This discusion shows that the graphs with two incommensurate 
lengths belong to a non-generic class because the singular manifold can never be 
of dimension n — 3 for rt = 2. It is important to notice that this statement does 
not contradict Berry's theory because he considers general Hamiltonian systems 
where the levels are given by an equation like /(fc, li, I2) — (if we call li and 
I2 the two parameters that enter in his theory) which allows the existence of a 
cone in the (fc, ^i, I2) space, while, for graphs with two incommensurate lengths, 
the secular equation has the special form f{k,li,l2) — f{kli,kl2) — which 
does not allow the existence of such a cone. 

In summary, the important point which makes the difference in the behavior 
of P(A) at small spacings A is the dimension of the singular manifold which 
is n — 3 for the repulsion and a distribution as P(A) ~ A, but n — 2 for the 
clustering and a distribution as P(A) ^ constant. 

We notice that the degree of the polynomial that describes the surface around 
the singular point seems to be not essential for this matter because we have 
repulsion and clustering for cases where the surface around the singular point 
is given by a quadratic form [see Subsection 5.1 Subsection 5.2 and Eq. (^)]. 
We have also seen an example with clustering [see Eq. (^] where the surface 
is described by a cubic form. 
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We have not commented on the topology of the singular manifold. It can 
happen that the intersection of two surfaces is transverse or tangent and this 
will influence the behavior of -P(A) near A = 0. In all the examples that we 
consider, the intersection is transverse which is the generic case when there is 
no restriction on the kind of surface. 

Several extensions of this work are possible. We may wonder which are the 
generic properties of the surface S and the function t(^) for typical graphs and 
apply the formula (|2^) to such generic situations. We notice that there are 
some restrictions on E. In particular, the number of sheets with projection in 
a given direction depends on the number of bonds with lengths associated with 
this direction, as we have seen in Section |[ There can be other restrictions 
as a consequence of the properties of the Hamiltonian (hermiticity, etc). In 
this problem, a difficulty comes from the fact that the dimension of the torus 
is equal to the number of incommensurate lengths in the graph. This makes 
difficult the study of complex graphs. Another possible direction is the study of 
perturbations of integrable graphs. Such perturbations are expected to deform 
the surface E and we may investigate the transition from our quasi-Poissonian 
distribution (|9|) to the RMT distribution. 

Since the graphs of the kind that we have discussed here have been used to 
model transport in mesoscopic systems, our work can find interesting applica- 
tions in this context. 



Appendix 

Here, we derive Eq. (^) from Eq. (^7|). The strategy is to divide the surfaces 
Efe into regions where the minimum, which appears in the definition (Uq) of 

Tk{sk), takes a given form. We call TZj the region where Tk{sk) = —f-^ with 

j ^ k and 7?.| the region where Tk{sk) = ^■ 
First, we compute 



/ dsi S[s - n{si)] = dxl--- dxl 5[s ~ Ti{xl, . 
Jt,i Jo Jo 

with Ti(si) = min^^i (ji^, 

Consider the fiow (Q) introduced in Section || 

xi = ht 

Xi = lit + x'^ , i = 2, . . . ,n (50) 

We look after any region in the surface Ei where ti(x2, . . . ,x^j) = j^. To 
determine these regions, we replace i by in (|5^). This region should satisfy 
the following inequalities 

< = /j-^ + < TT , i = 2,...,n (51) 
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which express the fact that the trajectory did not cross any boundary of the 
torus before arriving aX xi — tt. 

From we get that the region TZl where ti(si) = is given by 

0<x° <^(^1- , i^2,...,n (52) 
We note that this result gives a border for aU the other regions 

>7r(^l-|^^ , i = 2,...,n. (53) 

Now, we look for the regions where ti(si) = — j— ^, Vj ^ 1. Again, from the 
substitution of this expression in Eq. (pO|), we obtain the following inequalities 

< xi = ^-^{tt - x") < n 

Q < Xi = -^(tt — x") + < TT , z = 2,...,n, i^k 
'■J 

These inequalities imply that 

TT (^1 - 1^ < < TT , J ^ 1 

0< 4, <7r-^(^-a:°) , m^{l,j} (54) 

is the region TZj where ti(si) = ^ ' . 

The union of TZ\ given by (|2|) with TZj (Vj 7^ 1) given by ( p^ is equal to 

Si. 

Thus, we have that 
/ dsi 6[s - ri(si)] = 

The first term represents the integration over TZ\ and the second term the inte- 
gration over TZ^j. The explicit evaluation gives: 



n 



7r 1- 



dx' 



(55) 



dsi (5[s — ''■i(si)] = tt'' 



n 



n (- 



e 



(56) 
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Now, we compute 



dsk S[s - Tk{sk)] = dxl-- - j dxl_^ / dxl^^ ■■■ dx^ S[s - Tfc(a;?, . . . ,xl_^,xl^^, . 
Sfc Jo Jo Jo Jo 

(57) 

As before, we consider the flow 



Xk = ht 



, n, i ^ k 



(58) 



It is easy to see that there is no region where Tk{sk) — f^- This is due to the fact 
that the existence of such a region requires |^7r + x\ < it and, because li > Ik, 
there is no positive value of Xi where this inequality holds. 

Now, we look for the regions TZ^ where Tk{sk) — —f-^. Replacing this Tk{sk) 
in (p8) we get the following conditions 



< ifc 
Q <Xi 



(tt - x^) < tt 
■(tt — x") + < TT , i=l,...,n, i^k 



(59) 



These conditions are satisfied only in the region TZ^ where Tk{sk) — —j- 
which is defined by 



TT 1 - 



< X'' < TT 



0< 4. <^-^(^-4) 



and 

(60) 
(61) 



Again, the union of the regions Tih gives Sfc. Note that these regions are outside 



the border given by (53) as one expected. 
Thus, one has 



dsk 5[s - Tk{sk)] ^ 



n 



or more explicitly 

P n 

/ dsk S[s - rfe(sfe)] = y2 I J 



n ('^-M 



TT ~ X a 



(62) 



(63) 
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Finally, we substitute the results (p6|) and (|63|) in (Wlh and we get 



Pis) = 



(5 ( s - — 



^ n n 



fe=l j^^fc 



n ('^-^^^) 

(64) 



which gives Eq. (^9|) when written in terms of the variable A — -^^s. 
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